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We treat the propagation of nucleon in nuclear matter by evaluating the ensemble average of the 
two-point function of nucleon currents in the framework of the chiral effective field theory. We first 



derive the effective parameters of nucleon to one loop. The resulting formula for the effective mass 
was known previously and gives an absurd value at normal nuclear density. We then modify it 
following Weinberg's method for the two-nucleon system in the effective theory. Our results for the 
effective mass and the width of nucleon are compared with those in the literature. 

PACS numbers: 11.30.Rd, 12.38.Lg, 12.39.Fe, 24.85. +p 



I. INTRODUCTION 



Chiral perturbation theory is the effective theory of strong interactions, based only on the approximate chiral 
' symmetry of the QCD Lagrangian and its assumed spontaneous breaking. While this theory is eminently successful 
in dealing with pionic processes [l|, Q, an d also to a great extent, with those of pions and a single nucleon its 
direct application to processes involving two or more nucleons is problematic [H . As first pointed out by Weinberg [f| , 
i' the contribution of graphs for such processes grow with increasing number of loops, in gross violation to the power 
JL , counting rule. This failure of power counting is, of course, only to be expected, as precisely the sum over such (infinite 
number of) loop graphs would give rise to the bound states, such as the deuteron and the virtual bound state in the 
C"| , two nucleon system. 

A better analysis of graphs can be made in the old fashioned perturbation theory [5| . Here it is the loop diagrams 
. ! containing one or more pure-nucleon intermediate states that lead to divergence of the perturbation expansion, while 
those with intermediate states containing at least one pion obey the power counting rule. Accordingly Weinberg [j| 
proposes to construct first the effective potential from those connected graphs of the T-matrix, that do not contain 
any pure-nucleon intermediate state. Then the graphs with pure-nucleon intermediate states are summed over by the 
Lippmann-Schwinger intergal equation for the T- matrix with this effective potential. 

In this work we apply this procedure to find the effective mass of nucleon in nuclear matter. The one-loop formula 
is long known [6[, giving the mass-shift as the nucleon number density times a constant, depending on the parameters 
of the effective Lagrangian to leading order. When expressed in terms of NN scattering lengths, this constant and 
hence the mass-shift becomes too large to be acceptable. 

What goes wrong with the one-loop result is clearly the approximation of the appropriate NN scattering amplitude, 
that is present in the formula, by merely the s-wave scattering lengths. To improve this amplitude within the effective 
theory, the Weinberg analysis suggests that we regard the constant AW amplitude as the effective potential to leading 
order and solve the Lippmann-Schwinger equation for the complete amplitude. Replacing the constant amplitude in 
the one-loop formula by this solution, we get a greatly reduced value for the effective nucleon mass. We also find the 
effect of including phenomenologically the effective ranges in the formula Q • 

We first derive in some detail the old result for the nucleon effective mass, using the real time formulation of field 
theory in matter Q. We begin with the ensemble average of the two-point function of nucleon currents It is 
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evaluated to one-loop with vertices from the effective Lagrangian Q, using the method of external fields 0- The 
nucleon pole term gives immediately the the effective mass to leading order in the usual power counting. (We also get 
the residue at the pole, but it contains an unknown coupling constant.) The constant scattering amplitude contained 
in it is then improved upon as described above. 

In Sec. II we obtain the pieces of the effective Lagrangian, that will be needed in evaluating the one-loop graphs. 
We use the Dirac, rather than the Pauli, spinor for the nucleon. The advantage of this relativistic treatment is the 
clear separation of the vacuum and the density dependent parts of the nucleon propagator, the latter part containing 
the on-shell delta-function. In Sec. Ill we evaluate the graphs to obtain the leading terms for the effective parameters 
of nucleon in nuclear matter. We modify the mass shift formula in Sec. IV. Finally Sec.V includes some comments 
on our work. In Appendix A we describe a dispersion integral method to obtain the density dependent contributions 
from the one-loop graphs. The short Appendix B gives the result of partial wave expansion of the spin averaged AW 
scattering amplitude in the forward direction. 



II. CHIRAL LAGRANGIAN 

We study the nucleon propagation in nuclear matter by analysing the two-point function of the nucleon current 
v( x ) 01 j built out of three quark fields so as to have the quantum numbers of the nucleon. In this work we do not 
need to spell out the form of this current; it suffices only to note its tranformation under the chiral symmetry group 
SU(2) R x SU(2) L , 

VR * 9RVR-> Vl > 9lVLi g R eSU(2) R , g L e SU(2) L , (2.1) 

where the subscripts R and L on rj denote its right- and the left-handed components, t\r,l = |(1 ± 7s)^- 

To get the low energy structure of vertices with the nucleon current in the effective theory, we couple it to an 
external (spinor) field f(x), thereby extending the original QCD Lagrangian Cq CD to 

CQCD^C^cn+Jv + vf. (2.2) 

Writing in terms of R and L components, 

7v = 7b.Vl + 7lVR, 

we see that such a term is chirally invariant if the external fields f(x) R x transform oppositely to r](x) R ,L, 

Ir -> 9l/r, ft -» 9r/l ■ (2.3) 

The transformation rules for the pion and the nucleon fields are known 0, flf3 | . One can define three different quantities 
involving the pion triplet, namely U, u (it 2 = U) and = iu* (8^11)11' transforming as, 

U g R Ugl 

u — > gRiih} = hug L , 

hu^ , (2.4) 



u. 



where h is an element of the unbroken subgroup, h £ SU(2)y. The nucleon field ip(x) transforms according to the 
isospin i representation of SU(2)v, 

ip — ► hip. 

We can now construct the effective Lagrangian for the nucleon-nucleon system including pions in presence of the 
external field /. Its terms may be put into three groups, 

which we now write down one by one. is given by the familiar terms, 

= Zf-drUWrt + ^ - m)V> + \g A Hi^ • (2.6) 
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Here is the constant (<?a = 1-26) appearing in the neutron beta-decay. We choose the explicit representation 
U = exp(i<fi a T a j F n ), where <fi a are the hermitian pion fields (a = 1,2,3), r a , the Pauli matrices and F„, the pion 
decay constant (F^ — 93 MeV) . In the following we do not need vertices with pion fields only. 

The piece C^i giving the leading quartic interaction of nucleons has been written by Weinberg as Q , 

C r = -^-(N^Nf-^(N^Nf + --, 

where Cs and Ct are constants and N(x) is the Pauli spinor (and isospinor) for the nucleon and ' • • •' refers to terms 
with two and more derivatives. As we intend to work with the Dirac spinor ift(x), we rewrite it as 

V = W + 7oM 2 - ^pMl + 70)775^}' + • • • • (2.7) 
In the piece Cf dependent on the external field, we need two couplings, namely, that are linear and cubic in i/;, 

C f = Cfip + . 

Using the above transformation rules, the piece linear in ip may be written as 

c H = xJuJr)^ + avttjv + h.c. 

= 2/(1-75)^+ y7(l+7isW + ^c., (2.8) 
where invariance under parity requires A = A' . Following the construction of £.04 , we write the other piece as 

As - — At - — 

£fil>3 = —^(1 + 70)^/(1 + 70)^ + —^(1 + 7o)775^/(l + 70)775^ + h.c. (2.9) 

Two other possible non-derivative terms involving f do not appear here for the same reason as they did not in Eq.(2.6), 
namely that the anticommutativity of the nucleon field allows each of these terms to be written as linear combinations 
of the two terms written above. Clearly here also there are terms with two or more derivatives on the nucleon field. 



III. ONE-LOOP FORMULA 



We start with the ensemble average of the two-point function of nucleon currents in symmetric nuclear matter, 

n(?) =i [ d 4 xe lqx (Tr](x)fj(0)), (3.1) 



where for any operator O, (O) = Tr [pO] /Trp , p = e~P( H ~^ N \ H being the Hamiltonian of the system, (3 the 
inverse of temperature and N the nucleon number operator with chemical potential \i. In general, we deal here with 
functions of two variables, namely qo = E and \q\. 

In the real time version of quantum field theory in a medium, the two point function assumes the form of a 2 x 2 
matrix. But the dynamics is given essentially by a single analytic function, that is determined by the 11-component 
itself. Denoting the 11-component of the matrix amplitude by Fn(E, q), the corresponding analytic function F(E, q) 
has the spectral representation in E at fixed q Q , 

It cotoME'-rimirnF^j)^ 
ir J E' — E — iE'e 



where the integral runs over cuts of F\x(E, q). 




FIG. 1: Free propagator and one loop Feynman graphs for the two point function in the low energy region. The wavy line 
represents the nucleon current, while the continuous and the dashed internal lines are for nucleon and pion propagation. 



The simplest example of a two-point function in nuclear medium is, of course, the free nucleon propagator. Its 
11-component is given by (ri| . 1 



\s(p) n = I d*xe^{Tij,(x)il>(0)) n = W + m) 



2 

P 



{n- (uj p )e(po) + n + {uj p )6(-p )} 2n8(p 2 - m 2 ) 



(3.3) 



where are the distribution functions for nucleons and antinucleons, n^(ujp) = je^^pTM) + 1} , uj p = yj m 2 + p 
The corresponding analytic function S(p), as given by Eq. (3.2) is 

S(p) = -{p' + m) — 1 , 

jr — m z + te 

which is independent of and identical to the free propagator in vacuum. 

In the following we shall work in the limit of zero temperature, when the distribution functions become, n~ 
9(fi — u>p), n + — > 0. Then the nucleon number density, n is given by 

d 3 P af 2Pf 



,2 



r6([i-u} p ) 



(2tt) 3 ^ 37T 2 ' 

where pf is the Fermi momentum related to the chemical potential [i by /i 2 = m 2 + p 2 F . We want to calculate the 
density dependent part of H(q) in the low energy region to first order in n. To this end we draw all the Feynman 
graphs with one loop containing a nucleon line. In addition to the single nucleon line forming a loop, we include also 
the loop containing an additional pion line to account for singularities with the lowest threshold. They are depicted 
in Fig. 1 along with the free propagator graph (a). 



1 Going over to the non-relativistic limit and omitting anti-particle contribution, it reduces to the more familiar form Il2f . 

c / \ 6>(p - p F ) 6{p F - p) 
-Su(p) = — + , 

po — uj p + le po — LO p — it 
where the vacuum and the medium contributions are mixed up. 
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Before we work out the graphs, it is possible to simplify the piece CtM in the Lagrangian, needed for vertices in 
graphs (c) and (d) of Fig.l. As two of the nucleon fields are contracted at the vertex itself, we can do so already in 
Cf^s and retain only the density dependent part in the contraction. Thus referring to Eq.(2.2), we get from Eqs. (2.8 
- 9) the complete expression for r\ to leading order in the effective theory as, 

{27T ' T I 3 

l + ^75 + Cn(l+ 7 o)|^, ( = —(A s -A T ). (3.4) 

We shall work at q = 0, when the free propagation graph (a) takes the form 

U(E) (a) = X 2 S(E) = - X \. 1(1 + 7o ), (3.5) 
w hj — m + te 2 

in the vicinity of the nucleon pole. We now calculate the corrections to it produced by the rest of the graphs. 

To evaluate the self-energy diagram (b), we rewrite the four-nucleon interaction term given by Eq. (2.6) conveniently 

as 

V = -g E CP i AB T i CD (MB)$c1>D), (3.6) 

i=S,T 

where V s = (1 + 70) and T T = (1 + 70)775. Then the two point function for this diagram is given by 

11(E)(6) = -X 2 S(E)aS(E). (3.7) 

Here the self-energy a is a constant given by 



° = -\Y. C ^\ 7?T4 {-2MS(p) n r)r + rStfuT*}, (3.8) 

where the ir(ace) is over 7-matrices. On inserting the density dependent part of Su from Eq. (3.3), a can be 
immediately evaluated to give 

o \2 1 1 

IL(E) {b) = - — {C S ~ C T )n {E _ m)2 -{l + 70 ). (3.9) 

Note the absence of a simple pole in this contribution. From the constant vertex graphs (c) and (d), we get 

n(£0(o)+(«o = -4A 2 Cn-^l(l + 70). (3.10) 
E — m 2 

The graph (e) has the two-particle (tt and N) intermediate state. It is evaluated in detail in Appendix A. As 
expected, it does not give rise to a pole at E = m. Following similar steps, we get the contribution of the remaining 
graphs as 

3\ 2 g A n 1 1 
16mF2 E-m2 [ 



n (£)(/R( 9 ) - ^ o( 1 +7o) + ---, (3.11) 



, , 3\ 2 q 2 ,n 11,, , 

n(^) W - - J F \ {E -m)^ 1 + ^ + - (3 ' 12) 

where the ellipses denote non-pole terms. Again observe that the graph (h) does not give rise to a simple pole. 

Collecting the results for the simple and the double poles at E = m, we find the vacuum pole (3.5) to be modified 
in nuclear medium to 

A* 2 1 , 

o(l + 7o), 



E-m* +ie 2 

where 



AU-(-^-2 C ]ny ) 13-131 



8mF 2 



3 



The constant (Qs — Ct) can be related directly to experiment, but no such relation appears to exist for £. As is 
known already [6j and we shall also see below, the formula for m* is unacceptable at normal nuclear density. 
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FIG. 2: Nucleon-nucleon scattering amplitude in the tree approximation. 



IV. MODIFIED FORMULA 



If we note that the nucleoli mass-shift formula (3.14) is given by graphs (b) and (h) of Fig.l and recall the presence 
of the mass-shell delta-function in the density dependent part of the nucleon propagator (3.3), it is easy to guess that 
the coefficient of n in this formula must be related to some on-shcll NN scattering amplitude at threshold. To obtain 
the actual relation, we evaluate the scattering amplitude in the same chiral Lagrangian framework as that Sect. II. 
In the tree approximation the contributing Feynman graphs are shown in Fig. 2. 

We choose to calculate the spin averaged amplitude in the forward direction, 

M(pi,p 2 -> PuPz) = j ^2 M(p 1 ,a 1 ;p2,a 2 Pi,cr 1 ;p 2 ,a 2 ) ■ 

The propagating nucleon, say, a proton, may scatter with a proton or a neutron in the medium. Thus the isospin 
structure of M in a symmetric medium must be given by 

A4 M pp — >pp -)- Mpn — >pn ■ 

With this spin and isospin structure, the amplitude M as calculated from the graphs of Fig. 2 is [Q = p 2 — Pi), 
M= | C l {-trT i (p' 1 +m)T l (p' 2 +m) + 2trT l (p' 1 + m)-trT l ^ 2 + m)} 



i=S,T 

\ 2 

trfa + m)$y 5 (#2 + m)$~f 5 /(Q 2 - ml 



1 (JaA 
4 \2F*J 



i 

2 _2 ^ 



^ \ (c s - c T) — -;^-" v + (4.D 



_ (Ei + m){E 2 +m) g\ 
Am 2 

Ei and E 2 being the energies of the two nucleons. (The graph (b) in Fig. 2 does not contribute to M.) Thus at 
threshold M involves the same combination of constants as those in Eq.(3.14) for m* . 

As shown in Appendix B, the amplitude M at threshold can be expressed in terms of the s-wave spin-singlet and 
-triplet scattering lengths ai and to get 

37T 

to* — m = — (ai + 0,3)71. (4-2) 
2to 

With their experimental values as quoted in the Appendix and at normal nuclear density (pp = 1.36 /to -1 ), it gives 
to* — to = —620 MeV [6j, which is unacceptably large. 

To look for the source of the problem, it is suggestive to rewrite (3.14) and (4.2) as 



1 



n> - n, = -- I {-6m 2 (c s - C T + j%\ ) = -- [" ^ {-127rm(ai + a 3 )} (4.3) 

to J (27r) J 2m [ V 4F*J J to J (2tt)^2to 
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FIG. 3: Nucleon mass shift in nuclear matter as a function of Fermi momentum. The dashed and the solid curves are drawn 
with amplitudes f- ■ (k) and // (A) respectively. 

FIG. 4: Nucleon mass shift in nuclear matter as a function of Fermi momentum. 



At this point we recall Weinberg's method for the two-nucleon amplitude. The constants in the curly brackets 
in the integrands above represent the scattering amplitude in terms of the s-wave scattering lengths. Clearly this 
approximation is too bad in the low energy region, as it does not reproduce the nearby poles due to the bound and 
virtual states. But we can treat these constants as effective potentials (to lowest order) and solve the Lippmann- 
Schwinger equation with it to get a better representation of the partial wave amplitudes. As Weinberg himself shows 
H, the solution in this case turns out to be the unitarized version of the potential. Accordingly we replace the 
constants by the corresponding momentum dependent amplitudes satisfying (elastic) unitarity, 



fl 1) (k)= ( ---<*) , i= 1,3. (4.4) 



1 



a; 



Here k is the centre-of-mass momentum, to be distinguished from p, which is the momentum of the in-medium nucleons 
in the rest frame of the nucleon under consideration, the two being related by k 2 = m{^/ m 2 + p 2 — m)/2. 

An even better approximation of the scattering amplitude is the effective range approximation, that would result 
from including derivative terms in the effective theory 0] ■ The corresponding replacement would read 

- ai ^ fP(k)=(-±- + ±nk 2 -ikj , i=l,3, (4.5) 

where ri,3 are the effective ranges in the spin-singlet and -triplet s waves. Without trying to relate the effective ranges 
to the coupling constants of the higher order terms in the effective Lagrangian, we shall take their values as well as 
those of the scattering lengths from experiment. 

We thus get both the real and the imaginary parts of the pole position as, 

rn*-i, = m -^f^ { A(k) + M k)}- (4-6) 

where fi(k) stands for /! (k) or f^ 2 \k), according as we choose the replacement (4.4) or (4.5). The width 7 is to 
be interpreted as the damping rate of excitations with quantum numbers of the nucleon. The numerical evaluation 
of Eq.(4.6) is shown in Figs. 3 and 4. In normal nuclear matter, we get Am = m* — m = —33 MeV and 7/2 = 110 
MeV, with amplitudes in the effective range approximation (solid curves in Figs. 3 and 4). 

Our results may be compared with those in the literature. There have been a number of calculations of the nucleon 
spectral function in terms of the (off-shell) NN scattering amplitude evaluated with available NN potentials. We pick 
out from these works the real and the imaginary parts of the on-shell nucleon self-energy in normal nuclear matter at 
zero three momentum. Thus Baldo et al |13J give Am = —55 MeV, 7/2 = 63 MeV (as reproduced in [16(]). Benhar et 
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FIG. 5: Nucleon width in nuclear matter as a function of Fermi momentum. The dashed and the solid curves are drawn with 
the same amplitudes as in Fig. 3. 

al give Am = —68 MeV, 7/2 = 16 MeV. Jong et al give Am = —65 MeV, 7/2 = 25 MeV. Also a phenomenological 
determination of the width in terms of the differential cross-section for the AW scattering gives 7/2 > 25 MeV (lq |. 
Finally we compare with a calculation by us [l7| based on the virial formula for the nucleon pole position, which gives 
Am = -37 MeV, 7/2 = 112 MeV. 



The spectral function of nucleon in nuclear matter depends on the interaction in a two-nucleon system. To treat 
this interaction in the framework of effective chiral field theory [f| , one has to derive first the effective potential from 
the effective Lagrangian and then solve the dynamical equation with this potential. In this way one can restore the 
otherwise invalid power counting rule and accomodate at the same time the singularities of the scattering amplitude 
in the low energy region. 

The failure of the earlier calculation of the effective nucleon mass [(| may now be understood in this framework as 
due to representing the two-nucleon scattering amplitude by the (leading term of the) effective potential itself. In this 
work we proceed further to carry out the next dynamical step, replacing the potential by the unitarised scattering 
amplitude. By including the effective range terms in the s-wave amplitudes, we actually include the next-to-leading 
order terms also in the effective potential. 

It will be noted that our mass-shift and width formulae do not involve any Pauli blocking effect, our improved 
scattering amplitude being still in vacuum and not in medium. An estimate of this effect can be made from the work 
of Ref. [10|, where one finds a change of about 20% in the value of the potential energy per nucleon in nuclear matter, 
when the Pauli projection operators are withheld from the calculation. Incidentally, the phenomenological formula 
for the width in Ref.[l6| involves the differential cross-section in vacuum, the Pauli blocking factors appearing only 
for the final nucleons. Such factors, however, do not appear in our formula, as it depends on the forward scattering 
amplitude, where the final particles occupy the states vacated by the initial ones. 

The formula (4.6) bears a close resemblance to the one obtained from the virial expansion of the nucleon self-energy 
to first order [19j . In the virial formula, it is the full (spin averaged) forward amplitude that appears in the integral, 
which was evaluated in Ref. [13] using the phase shift analysis of experimental data on AW scattering [2(J, while the 
integral here includes only the part of this amplitude corresponding to s-waves in the effective range approximation. 
But in the low energy region over which the integrals arc evaluated, this approximation for the s-waves agrees well 
with the phase shift analysis. Further, the higher partial waves contribute negligibly to the virial formula. We thus 
expect the close agreement of the present results with those from the virial formula, stated at the end of the previous 
Section. 



V. 



DISCUSSION 
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appendix A 



Here we describe a method to evaluate the density dependent part of the amplitude for the one-loop graphs. In this 
method we calculate the imaginary part of the Feynman amplitude and construct its real part by a dispersion integral 
over it. Since we exclude the vacuum contribution and consider only the density dependent part, no subtraction 
is necessary for these integrals. The same result can be obtained in a simpler way by retaining only the density 
dependent part of the nuclcon propagator in the amplitude and using the mass-shell delta function in it to integrate 
out the energy component of the four-momentum. But the dispersion method is more transparent in that it shows 
the cut structures, where the contributions come from. 

We now evaluate graph (e) of Fig.l in some detail, it being the prototype of all other graphs. It has the amplitude 



where 



o\2 

n (9)(e) — ^r( 9 ) 



T(q) n =i I d^e^iflt^ns-SWuTs. 



(A.l) 



(A.2) 



Due to the presence of 0(±po) in Eq.(3.3) for S(p)n, it is not convenient to transform the amplitude in momentum 
space. Instead we integrate out the po variable in S(x)n itself, getting 



-S(x)n 

i 



d 3 p 



{(j* + m)(l - n-)e- lp - x + (jj - m)n+e lp - x } 9(x Q ) 

, Po = u p - 



(27r) 3 2w p 

{{$ + m)n-e- ip - x + (jj - m)(l - n + )e ip - x } 6{-x ) 



The analogous expression for the pion propagator is 



1 



tD{x)i 



d 3 k 



(27r) 3 2w fc 



{(1 + n)e- ik - x + ne ikx } 9(x ) + {ne~ lk - x + (1 + n)e ik - x ) 9(-x ) 



(A.3) 

fc = uj k (A.4) 
and n(u!k) are zero 



where n is the pion distribution function, n = (eP u}k — 1) — 1 , LUk — y rn^ + fc 2 . Though n + (ui. 
for the medium we are interested in, we retain them at this stage for generality and symmetry. 

With these expressions for the propagators, we may carry out both the x° and x integrations in Eq.(A.2) giving 
the energy denominator and the 3-momentum delta function respectively. The imaginary part may now be read off 
and put in the form, 



ImT{q) n = ir tanh(/3( 9o - n)/2)ImC(q), 



(A.5) 



where 



ImT{q) = J { J )3 P 2 ^ J { J )3 k ^ (jt> - m) { (1 - n + n)5^ (q - p - k) + (n~ + n)5^ (q-p+k)}. (A.6) 

Here the first term corresponds to rj — > ttN and the second to r/ + tt — * N, giving rise respectively to the discontinuity 
across the unitary and the 'short' cuts [2l[. With qo = E and for q — 0, at which we work, these cuts extend in the 
-B-plane over E > m + m v and < E < m — m OT . Note the opposite sign before n~ in the two terms. 
The complete result for ImT has also another piece given by 



d 3 p 



d 3 k 



(2tt) 3 2uj p J (2ir) 3 2uj k 



i (1 - n+ + n)8 {i) (q+p + k) + (n+ + n)8 (4) {q + p - k)\ , 



which is non-zero for E < only and has no term proportional to n~ . Clearly this observation does not depend on 
the structure of vertices in the loop diagrams. We thus have the general result that there is no spectral function to 
one loop for E < in a medium with nucleons only. 

The 3- momentum integrations in Eq.(A.6) can be carried out immediately to get 



ImT(E) = ± 



f(E) 
E ' 



f(E) = 



8tt 2 



-(70W -m)n (u>), 



(A.7) 
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where u> = {E 2 + m 2 — m 2 r )/2E is the nucleoli energy u> p , as restricted by the delta functions. The ± signs correspond 
to the unitary and the short cuts respectively. Inserting this result in Eq.(3.2), we get the spectral representation for 

r(£), 

/■°° dE'fjE') dE'f(E') 
T[E) - J m+m „ E'(E'-E) J E'(E>-E) ■ (A - 8) 

The range of the second integral can be mapped on to that of the first by the inverse transformation E' — > (to 2 — 
TO 2 )/i?'. Noting that u> and hence f(E) is form invariant under this transformation, we have 

Once we are on the unitary cut, the kinematics is determined by the first <5-function in Eq.(A.6). For q = 0, it gives 



E' = ^m 2 +p 2 + ^ml+p 2 , (A.10) 

so that E' and p are the total energy and the 3- momentum in the centre-of-mass frame of the intermediate n N system. 
The distribution function restricts the upper limit of the integral to the energy given by Eq.(A.lO) with p replaced 
by the Fermi momentum pp. Now setting = 0, Eq.(A.9) gives T to leading order as 

1 f PF dv v 2 

which is finite at E = m. 



appendix B 

The result (4.1) corresponds to s-wave amplitudes. To get the exact combination of these amplitudes, we recall the 
partial wave analysis of M. Following the usual convention, we write it in the center-of-mass system as 

M = S«w(p^ 1 \k) + l7^°\k ) y (B.l) 

where W and k are the total energy and momentum in the cm. frame. Each of the isospin amplitudes above has the 
partial wave expansion 22], 

7 (/) (fc) = 2.i^(2 J + l)/^(fc). (B.2) 

Here the factor 2 is due to the identity of the scattering particles. The total angular momentum j is obtained by 
coupling the total spin and orbital angular momenta s and /. In general, / /js is a 2 x 2 matrix in / space, whose 
diagonal elements enter the sum in Eq.(4.3). The quantum numbers of the partial wave amplitudes are governed by 
the antisymmetry of the total wave function, 

(-l)*.(-l)i-«.(-l)i-' = -l. 

Thus in the s-wave approximation, Eq.(4.2) reduces to 

M = 6TrW{h(k) + f 3 (k)} thr ^h old -127rm( ai + o 3 ), (B.3) 

where /i,3 are the spin singlet and triplet s-wave amplitudes and 01,3, the corresponding scattering lengths. Exper- 
imentally ai = —23.74 fm and 03 = 5.31 fm (23|. We also need the experimental values for the effective ranges, 
ri = 2.7 fm and r 3 = 1.70 fm ||. 
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